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ABSTRACT 

In the standard Cold Dark Matter (CDM) theory of structure formation, 
virialized minihalos (with T vir < 10 4 K) form in abundance at high redshift 
(z > 6), during the cosmic "dark ages." The hydrogen in these minihalos, the 
first nonlinear baryonic structures to form in the universe, is mostly neutral and 
sufficiently hot and dense to emit strongly at the 21-cm line. We calculate the 
emission from individual minihalos and the radiation background contributed 
by their combined effect. Minihalos create a "21-cm forest" of emission lines. 
We predict that the angular fluctuations in this 21-cm background should be 
detectable with the planned LOFAR and SKA radio arrays, thus providing a 
direct probe of structure formation during the "dark ages." Such a detection will 
serve to confirm the basic CDM paradigm while constraining the shape of the 
power-spectrum of primordial density fluctuations down to much smaller scales 
than have previously been constrained, the onset and duration of the reionization 
epoch, and the conditions which led to the first stars and quasars. We present 
results here for the currently-favored, flat ACDM model, for different tilts of the 
primordial power spectrum. 

Subject headings: cosmology: theory — diffuse radiation — intergalactic medium 
- large-scale structure of universe — galaxies: formation — radio lines: galaxies 



1. Introduction 

No direct observation of the universe during the period between the recombination 
epoch at redshift z ~ 10 3 and the reionization epoch at z > 6 has yet been reported. 
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While a number of suggestions for the future detection of the reionization epoch, itself, have 
been made, this period prior to the formation of the first stars and quasars - the cosmic 
"dark ages" (e.g. Rees 1999) - has been more elusive. Standard Big Bang cosmology in the 
CDM model predicts that nonlinear baryonic structure first emerges during this period, with 
virialized halos of dark and baryonic matter which span a range of masses from less than 
10 4 M Q to about 10 8 M Q which are filled with neutral hydrogen atoms. The atomic density 
n H and kinetic temperature T K of this gas are high enough that collisions populate the 
hyperfine levels of the ground state of these atoms in a ratio close to that of their statistical 
weights (3:1), with a spin temperature Ts that greatly exceeds the excitation temperature 
T* = 0.0681K. Since, as we shall show, T s > T CMB , the temperature of the Cosmic Microwave 
Background (CMB), as well, for the majority of the halos, these "minihalos" can be a 
detectable source of redshifted 21-cm line emission. The direct detection of minihalos at 
such high redshift would be an unprecedented measure of the density fluctuations in the 
baryons and of the total matter power spectrum at small scales, which will not be probed 
by other methods yet discussed (e.g. CMB anisotropy). 

The possibility of 21-cm line emission or absorption by neutral H at high redshift has 
been considered before (Hogan & Rees 1979; Scott & Rees 1990; Subramanian & Padman- 
abhan 1993; Kumar, Padmanabhan, & Subramanian 1995; Bagla, Nath, & Padmanabhan 
1997; Madau, Meiksen & Rees 1997; Shaver, et al. 1999; Tozzi et al. 2000). Prior to the 
release of radiation by nonlinear baryonic structures which condense out of the background 
universe, the spin temperature of HI in the diffuse, uncollapsed gas in the intergalactic 
medium (IGM) is coupled to the CMB, so that T s = T CMB and neither emission nor ab- 
sorption in the 21-cm line is possible. Recently, attention has focused on the possibility that 
radiation from early stars and quasars might decouple Ts from Tcmb by "Lya pumping" - 
resonant scattering in the H Lja transition followed by decay of the upper state n = 2 to 
the ground state n = 1 into one or the other of the hyperfine levels (Madau et al. 1997; 
Tozzi et al. 2000). This mechanism, it has been suggested, will operate on HI in the diffuse, 
uncollapsed IGM during reionization, first to make T s < T CM b, so that the 21-cm transition 
can be seen in absorption against the CMB, until the same Lya scattering heats the gas 
shortly thereafter and makes T s > T C mb, thereby causing 21-cm emission in excess of the 
CMB, before reionization finally destroys the HI. In what follows, however, we show that 
a substantial fraction of the baryons in the universe may already have condensed out of the 
diffuse IGM into virialized minihalos, prior to and during reionization. Under these con- 
ditions, collisional excitation alone is sufficient to decouple Ts from Tcmb and cause 21-cm 
emission in excess of the CMB, thereby providing a signature of the cosmic "dark ages" and 
of their retreat during reionization. 
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2. 21-cm Emission from Individual Minihalos 

The 21-cm emission from a single halo depends upon its internal atomic density, tem- 
perature, and velocity structure. We model each CDM minihalo here as a nonsingular, 
truncated isothermal sphere ("TIS") of dark matter and baryons in virial and hydrostatic 
equilibrium, in good agreement with the results of gas and N-body simulations from realistic 
initial conditions (Shapiro, Iliev & Raga 1999; Iliev & Shapiro 2001, 2002). This model 
uniquely specifies the internal structure of each halo (e.g. total and core sizes, r t and r , 
central total mass density po, dark matter velocity dispersion ay = (47rpo r o) 1,/2 ; an d g as tem- 
perature T K = fim p <T V /kB, where \x is the mean molecular weight), for a given background 
cosmology as functions of two parameters, the total mass M and collapse redshift z co \\. 

The minihalos which contribute significantly to the 21-cm emission span a mass range 
from M m i n to M max which varies with redshift. M m \ n is close to the Jeans mass of the 
uncollapsed IGM prior to reionization, Mj = 5.7 x 10 3 (fi ^ 2 /0.15) _1/2 (Q b h 2 /0.02)~ 3/5 
x [(1 + ^)/10] 3/2 M Q , while M max = 3.95 x 10 7 

x (fi /i 2 /0.15)- 1/2 [(l + z)/10}~ 3/2 is the mass for which T vir = 10 4 K according to the TIS 
model (Iliev & Shapiro 2001) (since halos with T vir > 10 4 K are largely collisionally ionized). 
Halos with T vir > 10 4 K may have radiatively cooled gas inside them which would add to 
the signal we compute, but such gas is expected to lead to the formation of internal sources 
of ionizing radiation which will largely offset the effect. Since these additional effects are 
highly uncertain and are related to the onset of radiative feedback and reionization which 
we are neglecting in these calculations, we will not consider the role of higher temperature 
halos further. 

The flux per unit frequency, $ u = (dF / dv) reci received at redshift z = at frequency 
z/ rec from a minihalo at redshift z which emits at frequency u cm = v rec (l + z) is expressed in 
terms of the brightness temperature T^em = T^ rec {l + z) according to 

2z/ r 2 cc 

3Wec — p^B^VeclAf^halo , (1) 

where (Af2) halo = nrf/D\ = n (A9 halo /2) 2 is the solid angle subtended by the minihalo, and 
Da is the angular diameter distance. The brightness temperature T feem is determined by 
solving the equation of radiative transfer to derive the brightness profile of the minihalo and 
integrating this profile over the projected surface follows. 

The brightness temperature along a line of sight thru a minihalo at projected distance 
r from the center obeys the equation 

T b (r) = T CMB e" r(r) + f ^ Tge^'dr' , (2) 
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where quantities are defined in the comoving frame of the minihalo, frequency v refers 
here and henceforth to u cm , r(r) is the total optical depth thru the halo, and the effective 
absorption coefficient k u is given when T* <C T s by: 

3c 2 A l0 n m ,, N T* 
k v — 7, — ]\v — 3 

(Field 1958), where A w = 2.85 x 10 _15 s _1 is the Einstein A-coefficient for the 21-cm 
transition and /(z/) is the normalized line profile. The spin temperature Ts is determined 
by the balance between collisional and radiative excitation and de-excitation by atoms and 
electrons and by CMB and Lya photons, respectively, according to 

rp T C mb + VcTr + y a T a , . \ 

1 S = r- : , [V 

l+Vc + Va 

where T a is the color temperature of the Lya photons, and y a and y c are radiative and 
collisional excitation efficiencies, respectively (Purcell & Field 1956; Field 1958, 1959). The 
efficiency y c includes contributions from H° — H° collisions, y# , and from e~ — H° collisions, 
y e . Prior to the reionization epoch, Lya pumping is unimportant and collisional excitation 
alone must compete with excitation by the CMB. This is only possible for gas that is highly 
nonlinear and sufficiently hot. Such conditions are achieved only inside virialized halos. As 
shown in Figure 1, the optical depth of an individual halo is not negligible, particularly 
for smaller-mass halos (due to their lower Ts). Since Ts varies with radial position inside 
the halo, as a result of its significant central concentration, we must integrate equation (2) 
numerically. The face-averaged TJ, of this single halo is given by (Tfe)haio = (/ Tb(r)dA)/A, 
where A(M, z) is the geometric cross-section of a halo of mass M and collapse redshift z. 
The observed flux from an individual halo is then expressed with respect to the CMB by the 
differential antenna temperature 5T b = [(T 6 ) halo — T CMB (z)]/(l + z). 

The line-integrated flux F(M, z) received from this minihalo is equal to the flux calcu- 
lated for v — vq multiplied by a redshifted effective line-width Av e f{(z), defined by Au e g(z) = 
(J dd^)/du - For an optically thin minihalo, Au eS reduces to Au cS (z) = [f(y )(l + z)]^ 1 . In 
that case, for a thermal-Doppler-broadened line profile, Au cS (z) = [(27r / u) 1 / 2 z/ o-y/c](l + 2;) _1 . 
We have checked that this approximation is adequate even for the optically thicker halos at 
the small-mass end of the mass function. The differential line-integrated flux 5F(M, z) is 
given by replacing T^ rec in equation (1) by 5Tb and integrating over frequency as described 
above. 

Our results for individual minihalos are summarized in Figure 1. Line profiles of different 
minihalos along the same line of sight should not typically overlap. The proper mean free 
path A m f p = (^haioOhaio) -1 for photons to encounter minihalos in ACDM is 160 kpc at z = 9 
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(Shapiro 2001), corresponding to a frequency separation, Az/ sep ps voH(z)X m f p /[c(l + z)\ ~ 
0.1 MHz > Az/ Cff < 10 kHz. 

These results predict a "21-cm forest" of minihalo emission lines. At z = 9, for example, 
there are about 160 minihalo lines per unit redshift along a typical line of sight in an untilted 
ACDM universe (Shapiro 2001). Detecting the stronger lines would require sub-arcsecond 
spatial resolution, ~ 1 kHz frequency resolution, and ~ nJy sensitivity. SKA is expected to 
have sufficient resolution for such observation, but probably not sufficient sensitivity. 



The average differential flux per unit frequency relative to that of the CMB from all the 
minihalos observed within a given beam of angular size A0 beam and frequency bin Az/ obs is: 



where d 2 V(z)/dzdfl is the comoving volume per unit redshift per unit solid angle, the solid 
angle (Af2) bcam = 7r(A# bca m/2) 2 , and Az/ obs /A,2 = vq/(1 + z) 2 . We calculate the comoving 
density of halos at different redshifts using the Press-Schechter (PS) approximation for the 
halo mass function dn/dM. If we define the beam-averaged "effective" differential antenna 
temperature 5T b using 5$ u = 2z/ 2 /c B 5T f) (Af2) bcam /c 2 , then 



We consider the currently-favored, flat CDM model with cosmological constant ("ACDM", 
tt = 0.3, A = 0.7, COBE-normalized, Vt b h 2 = 0.02, h = 0.7), for three values of the 
primordial power spectrum index n p = 0.9, 1, and 1.1, using the primordial power spectrum 
of Eisenstein & Hu (1999). 

Results for 8$ v and 5T b are plotted in Figure 2. In principle, the variation of 5T b with 
observed frequency implied by the redshift variations in Figure 2 should permit a discrim- 
ination between the 21-cm emission from minihalos and the CMB and other backgrounds, 
due to their very different frequency dependences. However, the average differential bright- 
ness temperature of this minihalo background is very low and its evolution is fairly smooth, 
so such measurement may be difficult in practice with currently planned instruments like 
LOFAR and SKA. The angular fluctuations in this emission, on the other hand, should be 
much easier to detect, as discussed in the next section. 



3. 21-cm Radiation Background from Minihalos 
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4. Angular Fluctuations in the 21-cm Emission Background 

The amplitude of q-a angular fluctuations (i.e. q times the rms value) in the differential 
antenna temperature is given in the linear regime by 

= Qb(z)a p , (7) 

0±b 

where a p is the rms mass fluctuation at redshift z in a randomly placed cylinder which 
corresponds to the observational volume defined by the detector angular beam size, A^bcam, 
and frequency bandwidth, Az/ Q b s , and b(z) is the bias factor which accounts for the clustering 
of rare density peaks relative to the mass. We assume b(z) is the flux-weighted average over 
the mass function of b(M, z) — 1 + [y\ — 1)/S C , the linear bias factor, where = 8 c /a(M), 8 C 
is the value of the linearly extrapolated value of overdensity Sp/p corresponding to the epoch 
when a top-hat collapse reaches infinite density, and cr(M) is the the standard deviation of 
the density contrast filtered on mass scale M (e.g. Mo & White 1996). For a cylinder of 
comoving radius R = A# b cam(l + z)Da(z)/2, and length + z)cH (z)^ (Av / v) ohs , we 

have : 

2 8D' 2 (z) f°° f 1 sm 2 (kLx/2)J 2 [kR(l-x 2 y/ 2 } 



c n = — „„ , / dk dx 



p n 2 R 2 L 2 J J x 2 (l-x 2 ) 

x(l + /- 2 ) 2 ^ (8) 

(Tozzi et al. 2000) (with several typos in the corresponding expression in that paper cor- 
rected here), where D(z) = 5 + (0)/5+(z) is the linear growth factor, P(k) is the linear power 
spectrum at z — 0, and the factor (1 + fx 2 ) 2 , with / ps [r2(z)] - 6 , is the correction to the 
cylinder length for the departure from Hubble expansion due to peculiar velocities (Kaiser 
1987). 

Illustrative results are plotted for 3-cx fluctuations as a function of A#beam, for z = 7 
and 8.5, in Figure 3, along with the expected sensitivity limits for the planned LOFAR 
(300 m filled aperture) and SKA (1 km filled aperture) arrays. We plot in Figure 4 the 
predicted spectral variation of these fluctuations vs. redshift z for illustrative beam sizes of 
A#beam = 9' and 25'. These 3-cr fluctuations should be observable with both LOFAR and 
SKA with integration times of between 100 and 1000 hours. For a 25' beam, for example, 
3-cr fluctuations can be detected for untilted ACDM by both with a 100 h integration for 
z ~ 6 — 7.5 and a 1000 h integration for z < 11.5, while for a 9' beam, SKA can detect 
them after 100 h for z < 9 and after 1000 h for z < 13. Results for different values of z and 
A#beam are available upon request. 



-7- 



Acknowledgments 

We are grateful to E. Scannapieco for clarifying the effects of bias and referee P. Tozzi for 
his thoughtful comments. This work was supported by European Community RTN contract 
HPRN-CT2000-00126 RG29185 and grants NASA ATP NAG5-10825 and NAG5-10826 and 
Texas Advanced Research Program 3658-0624-1999. 

REFERENCES 

Bagla, J. S., Nath, B., & Padmanabhan, T. 1997, MNRAS, 289, 671 

Field, G. B. 1958, Proc. I.R.E., 46, 240 

Field, G. B. 1959, ApJ, 129, 536 

Eisenstein, D. J., & Hu, W. 1999, ApJ, 511, 5 

Hogan, C. J., & Rees, M. J. 1979, MNRAS, 188, 791 

Iliev, I. T., & Shapiro, P. R. 2001, MNRAS, 325, 468 

Iliev, I. T., & Shapiro, P. R. 2002, in "The Mass of Galaxies at Low and High Redshift" (ESO 
Astrophysics Symposia), eds. R. Bender & A. Renzini (Heidelberg: Springer- Verlag) , 
in press (astro-ph/0112427) 

Kaiser, N. 1987, MNRAS, 227, 1 

Kumar, A., Padmanabhan, T., & Subramanian, K. 1995, MNRAS, 272, 544 

Madau, P., Meiksen, A., & Rees, M. J. 1997, ApJ, 475, 429 

Mo, H., & White, S. D. M. 1996, MNRAS, 282, 347 

Purcell, E. M., & Field, G. B. 1956, ApJ, 124, 542 

Rees, M. J. 1999, Physics Reports, 333, 203 

Scott, D., & Rees, M. J. 1990, MNRAS, 247, 510 

Shapiro, P. R. 2001, in Proceedings of the 20th Texas Symposium on Relativistic Astro- 
physics and Cosmology, eds. H. Martel and J. C. Wheeler, (AIP Conference Series), 
pp. 219-232 



- 8- 



Shapiro, P. R., Iliev, I. T., & Raga, A. C. 1999, MNRAS, 307, 203 

Shaver, P. A., Windhorst, R. A., Madau, P., & de Bruyn, A. G. 1999, A&A, 345, 390 

Subramanian, K., & Padmanabhan, T. 1993, MNRAS, 265, 101 

Tozzi, P., Madau, P., Meiksen, A., & Rees, M. J. 2000, ApJ, 528, 597 

Fig. 1. — Individual minihalo sources of redshifted 21-cm emission in ACDM, redshifts 1+z = 
7 (short-dashed line), 10 (solid line), 15 (long-dashed line), and 20 (dotted line) vs. total 
mass of minihalo M. From top to bottom: optical depth r UQ (r = 0) at line-centered frequency 
u thru minihalo center, differential antenna temperature 5Tb, line-integrated differential flux 
8F{M 1 z) relative to the CMB, total differential flux per unit frequency fo, , angular size of 
minihalo (A6>) ha i , and redshifted effective width Au cS (z) of the 21-cm line as observed at 
z = at received frequency z/ rec = z/ (l + z)~ x . 

Fig. 2. — Minihalo radiation background. Average observed differential antenna temperature 
5T b and average differential flux per unit frequency 5$ u for beam size of A# bcam = 10' at 
the redshifted 21-cm line frequency due to minihalos vs. redshift z for ACDM models with 
power-spectrum tilts n p = 0.9, 1.0, and 1.1, as labelled. 

Fig. 3. — Predicted 3-cr differential antenna temperature fluctuations at z = 7 ((z/ rcc = 177.5 
MHz; top panel) and z = 8.5 (z/ rcc = 150 MHz; bottom panel) for bandwidth Az/ obs = 1 MHz 
vs. angular scale A6* beam for ACDM models with tilt n p = 0.9, 1.0, and 1.1, as labelled (solid 
curves) . Also indicated is the predicted sensitivity of LOFAR and SKA for a confidence level 
of 5 times the noise level after integration times of 100 h (dashed lines) and 1000 h (dotted 
lines), with compact subaperture (horizontal lines) and extended configuration needed to 
achieve higher resolution (diagonal lines) (see http://www.lofar.org/science). 

Fig. 4. — Predicted 3-a differential antenna temperature fluctuations at A6> bcam = 9' (top 
panel) and 25' (bottom panel) vs. redshift z for ACDM models with tilt n p = 0.9, 1.0, 
and 1.1, as labelled (solid curves). As in Figure 3, we also plot the predicted sensitivity for 
integration times 100 h (dashed) and 1000 h (dotted) of both LOFAR ("L") and SKA ("S"), 
as labelled (for bottom panel, sensitivity curves for LOFAR and SKA are identical), for com- 
pact subaperture and assuming rms sensitivity oc v~ 2A (see http://www.lofar.org/science). 
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